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ABSTRACT 

Elementary results are presented for the buckling of stitched, laminated composites 
containing delamination cracks. The stitching fibers are assumed to provide continuous, linear 
restoring tractions opposing the deflection of the delaminated layer adjacent to the crack. It is 
shown that there exists a characteristic length a for buckling: if the length, 2a, of the 
delamination crack exceeds 2a Q , then, when buckling occurs, it will consist of waves of period 
2a 0 and will usually not span the whole delamination. Simple expressions are derived for the 
critical buckling load and the minimum stitching density required to suppress buckling of the 
delaminated layer. 
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INTRODUCTION 

One of the principal obstacles to using relatively cheap graphite epoxy laminates in 
the commercial aircraft industry is their susceptibility to delamination, especially during 
impact, and the subsequent catastrophic growth of the delamination crack when the 
delaminated layer buckles under in-plane compression [see, for example, ref. 1], One promising 
solution to this problem is the incorporation of fiber tows normal to the laminate plane by 
stitching. While stitching tows do not eliminate delamination during impact, they do minimize 
loss of strength under subsequent compression [refs. 2,3]. It appears that the stitching tows 
bridge the delamination crack and prevent or reduce buckling of the adjacent delaminated 
layer. The in-plane stiffness then survives relatively unimpaired and, if the delaminated layer 
remains flat, the delamination crack does not grow, since it experiences no driving force. 

Failure under compressive loading occurs by some other mechanism and the compressive 
strength after impact is restored to an acceptably high value. 

This paper presents the simplest possible description of the mechanics of this buckling 
problem, which can be modeled as that of a buckling plate on a Winkler foundation of damping 
springs [e.g., refs. 4,5]. However, in contrast to the usual assumption, buckling deflections in 
the present problem can occur in one direction only. This constraint gives rise to a characteris- 
tic length for buckling and an enhanced value of the critical force for buckling. Simple 
arguments are presented to show how these fundamental quantities vary with stitching fiber 
density and the thicknesses of the delaminated layer and the substrate beneath it. 

BUCKLING IN THE PRESENCE OF LINEAR DAMPING 

Consider a delamination crack of length 2a lying along the x-axis, as shown in 
Figure 1. Let the crack be in a state of plane strain. Suppose the body containing the crack is 
subject to a compressive load that results in the compressive force F per unit length being 
imposed on the thin layer of delaminated material above the delamination crack. Buckling of 
the thin layer will result in the deflection w(x), which can be determined if the layer is thin 
enough by consideration of the balance of forces according to the elementary theory of bending 
plates [refs. 5,6]. If the delamination crack is bridged by stitching fibers, as shown 
schematically in Figure 1, the problem is modified by the lateral tractions those fibers impose. 
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Fig. 1. Schematic of a delaminated and buckled layer, showing lateral tractions opposing the 
buckling deflection. 

If the stitching fiber spacing is appreciably smaller than 2a, the lateral tractions can be 
considered continuous over x. Furthermore, if the stitching fibers are debonded from the 
matrix or are much stiffer than the matrix, then the tractions, denoted q, will be a linear 
function of the deflection: 

q(x) = ew(x) , (B>0) . (1) 


The spring constant 0 will be related below to the properties of the stitching tows. The 
deflection profile is given by the linear differential equation [refs. 4-6] 



D dx^ D 


0 


( 2 ) 


where F > 0 denotes a compressive force and D is the flexural rigidity of the delaminated 
layer. (The body has been assumed to be isotropic and homogeneous, a crude representation of a 
stitched laminate but a useful simplification when exploring fundamental aspects of the 
problem.) The deflection must also satisfy clamped end boundary conditions and be positive, 
since negative deflections would imply interpenetration. Equation (2) has both symmetric and 
antisymmetric solutions, but only the symmetric solutions can satisfy w > 0. The symmetric 
solutions have the form 


w(x) = A + coss + x + A coss x 
where s + = Jf (1± \/l 


b/f 2 ) 


(3a) 

(3b) 


with f = F/2D and b s b/D 


(3c) 


The boundary conditions determine both the ratio A + /A_ and the critical buckling load f (a) 
for the m th buckling mode for delamination crack length 2a. Numerically determined values of 


:(m) 


fj-' (a) are shown for the first few modes in Fig. 2 (see also [Ref. 4]). Each is a 

monotonically decreasing function of a, in contrast to the case of a plate with simply supported 
ends, for which the buckling load decreases for small a and then increases as a 2 for large a [e.g., 
ref. 5]. 
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Fig. 2 The critical force for buckling for the first few symmetrical buckling modes as a 
function of delamination crack length. 


With f 
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Fig 

it contains m maxi 
through zero and w 


determined, the corresponding deflection w^ m \x) is also determined to 

a l/ift tKci r\r\inte ^ jj-j 


hin a multiplicative constant by the boundary conditions. To the left of the points 
. 2, w< m >(x) is of one sign, i.e. it corLtains no zeroes in (-a,a); and, if it is taken to be positive, 
™ At the points A', B' , etc., d 2 w' rn '/dx 2 evaluated at x = ±a passes 

(x) acquires two zeroes, which begin at x = ±a and move into the interval 


(-a, a) as a increases. At all values of a to the right of A^ m \ w^ m \x) possesses at least two 
zeroes, i.e. it is no longer of one sign. 

For w > 0, the minimum buckling load, f Q , is thus found at the point i.e. 


f 


o 


3 


The crack length corresponding to point 



is given by 


(4) 


✓ 3 TT 

0 " 2 b* 


(5) 


It denotes a characteristic buckling length : whenever a > a D , buckling will occur not over the 
whole delamination (-a, a), but over some subinterval or series of subintervals of length 2a Q . The 
buckling profile over each subinterval is that for E; + a = 3n/2 and 5_a = tt/2, i.e. 


w ( x ) = A[cos + 3cos — ] 

0 2a 2a 


( 6 ) 
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with A an undetermined constant. The appearance of buckling when a>a Q is therefore 
qualitatively different from the case of an unbridged delamination crack. If a<a D , buckling 
occurs over the whole interval (-a, a) and the critical load rises with decreasing a as in Fig. 2. 

STITCHED LAMINATES 


For stitched laminates, the stiffness parameter e depends on how the delaminated 
layer is coupled to the rest of the material. Figure 3(a) shows the case of a delamination crack 
lying in the mid-plane of a thin panel. Buckling upon compressive loading occurs in a symmetric 
manner. Figure 3(b) shows the case of a thin delaminated layer lying over a thick substrate, 
which remains straight while the delaminated layer buckles. If the stitching fiber tows pass 
from top to bottom of the panel and they are initially unstrained, then the stiffness parameter 
is approximately 


8 = 



(case I) 


(7a) 

(case II) 

9 

(7b) 


where case I refers to the case of Fig. 3(a) and case II to that of Fig. 3(b), u $ is the volume 
fraction of the stitching fibers, Ef is the fiber modulus, and the dimensions n and t are defined 
in Figure 3. 
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Fig. 3 (a) Symmetric buckling, (b) Buckling of a thin layer over a thick substrate. 

Now the results of Section 2 were obtained for an isotropic material, whereas a 
stitched laminate is more closely orthotropic. Nevertheless, useful results in terms of orders of 
magnitude can be obtained by making some crude approximations concerning elastic 
properties. Since buckling of thin plates is determined mainly by in-plane elastic properties, the 
flexural rigidity of Eq. (3) can be approximated by 
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( 8 ) 


D = 


“f E f h3 

12(1 - v 2 ) 


where v is Poisson's ratio, ut is the volume fraction of fibers lying in the axial direction (parallel 
to the x-axis and the applied compressive load), and it has been assumed that these fibers have 
the same modulus as the stitching fibers. The stitching and axial fibers are indeed often of the 
same kind. The characteristic buckling length can now be approximated by 

h (case I) 

s 

c AVll* h 

c l t -u s J l h J n (case II) 

3^ Tt 

with c, = Vx - 1.5 when v = 0.3 ; 

1 2 3/Z (1-v )^ 

and the critical force f 0 of Eq. (U) is approximately 

^(ufUg)* 5 E f h (case I) (l° a ) 

c 2< u f u s> % $>* E f h (case II) (10b) 




(9a) 

(9b) 

(9c) 


with 


C 2 " 


3 3 / 2 ( 1 - v 2 )* 5 


1.0 when v = 0.3 


(10c) 


Equations (9) and (10) illustrate in a simple way how geometry and fiber density 
control the efficacy of stitching in suppressing buckling. In particular, Eq. (9) shows that the 
characteristic buckling length in stitched materials will usually be no more than an_ order of 
magnitude greater than the thickness, h, of the delaminated layer. Even if u $ ~ 10” 2 Uf, i.e., 
sparse stitching, one still l^s a Q <* 5h for case I, since the ratio Uf/u s appears to the power 
The additional factor (t/h)" in case II is also unlikely to be much greater than 2 in practice. 

In the absence of impact damage, stitched laminates generally fail under in-plane 
compression by the formation of a kink band of buckled and broken axial fibers. The stress at 
which this occurs corresponds to a particular value, f ' , of the force acting on the ends of a 
delaminated layer if it has not yet buckled. Substituting f' for f Q , Eq. (10) provides a simple 
estimate of the stitching fiber density required to suppress Duckling and eliminate delamination 
crack growth as a potential failure mechanism. Thus the required value of is proportional 
t° (f^/h) 2 . Conversely, the critical force, f Q , is more sensitive to fiber volume fractions than 
is the characteristic buckling length, being proportional to u $ 2 . 

Equation (10) also shows that f Q in case II varies as t 2 , which simply reflects the fact 
that stitching fibers of shorter initial through-thickness length experience greater strain for a 
given buckling deflection. Thus an effective method of raising f Q for a given density of 
stitching fibers and laminate thickness is to pass stitching tows only part way through the 
laminate, achieving through-thickness reinforcement by staggering the stitching at different 
depths. Laminates with such stitching patterns are indeed available, but no results concerning 
their delamination and buckling behavior have yet been published. 
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CONCLUSIONS 


The presence of stitching tows introduces a minimum compressive load, f Q , required 
to buckle a delaminated layer in a stitched laminate, regardless of the length of the 
delamination crack; and a characteristic buckling length, a Q . If the delamination crack length 
2a exceeds 2a Q , buckling has the form of waves of length 2a Q , which arise when the load 
exceeds f^. In typical stitched laminates, the length a Q will not exceed the thickness, h, of the 
buckling delaminated layer by more than one order of magnitude. The critical force, f , can be 
enhanced by increasing stitching density or passing stitching tows only part way through the 
laminate. 
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